A Bose-Einstein condensate in a random potential 
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An optical speckle potential is used to investigate the static and dynamic properties of a Bose- 
Einstein condensate in the presence of disorder. For strong disorder the condensate is localized in 
the deep wells of the potential. With smaller levels of disorder, stripes are observed in the expanded 
density profile and strong damping of dipole and quadrupole oscillations is seen. Uncorrelated 
frequency shifts of the two modes are measured for a weak disorder and are explained using a 
sum-rules approach and by the numerical solution of the Gross-Pitaevskii equation. 
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The nature of a bosonic system in the presence of dis- 
order has been intensely explored both theoretically and 
experimentally, particularly in the context of ^He in a 
porous material, where the breakdown of superfluidity 
with sufficient disorder was observed 1]. In the zero 
temperature limit, theoretical studies have provided the 
striking result that there exists a regime where the su- 
perfluid fraction is significantly smaller than the conden- 
sate fraction 2]. A further consequence of disorder is 
the quantum phase transition from a superfluid to a lo- 
calized Bose-Glass state js^. Indications of such a state 
have been found in several physical systems but the pre- 
cise phase diagram of the progression from superfluid, to 
Bose-Glass, to Mott-insulator state remains unclear ^J. 
Furthermore, Anderson localization occurs for weakly in- 
teracting bosons 5] , and plays an important role in solid 
state physics where electron transport can be disrupted 
by defects in a solid , and more recently was observed 
for photons in strongly scattering semiconductor powders 

Q. 

The realization of Bose-Einstein condensation (BEG) 
in a dilute alkali vapor offers the possibility of an ideal 
system in which to explore the rich arena of disorder re- 
lated phenomena. Already there has been considerable 
work on the superfluidity and long-range coherence prop- 
erties of these degenerate Bose gases, and with the com- 
bination of BEG and optical lattices there has been an 
explosion of observations on phenomena of solid state- 
physics 0, notably including the superfluid to Mott- 
Insulator quantum phase transition y\. Recently, theo- 
retical works have discussed the possibility of Bose-Glass 
and Anderson localization transition in a BEG in a dis- 
ordered optical lattice 0|- Understanding the ef- 
fect of disorder is also important for BEG in microtraps 
where fragmentation of the trapped BEG density pro- 
file has been observed and is now attributed to intrin- 
sic disorder in the fabrication of the microchip 0, . 
It is crucial for integrated atom optics on microchips to 
quantitatively define how this disorder will change the 
coherence and transport properties, particularly consid- 
ering that in ID disordered systems localization of the 
excitations could play an important role 

We present in this Letter the first results on a Rb BEG 
in a disordered potential, where the precise realization of 
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FIG. 1: (a) Optical setup for both the speckle potential and 
the imaging beam for the BEC. The axial direction of the 
magnetic trap is in the vertical direction of the figure, (b) 3D 
representation of the speckle potential (left) and its Fourier 
transform (right). The dotted lines correspond to a length 
scale of about 10 fim in the axial direction. 



the disorder is controllable. We consider both the limit of 
strong disorder, where we observe that the BEG is local- 
ized in the wells of the high random potential, and the 
situation where the disorder is smaller than the chem- 
ical potential of the BEG, where we see stripes in the 
expanded density profile and damping of dipole oscilla- 
tions. In the perturbative limit of very weak disorder we 
observe shifts in the frequency of the quadrupole excita- 
tion mode. 

To produce the random potential a laser beam is shone 
through a diffusive plate. The beam is derived from a 
Ti:Sa laser at A = 822 nm, far detuned with respect to 
the Rb DI line at A = 795 nm. The resulting speckle 
pattern is imaged onto the trapped BEG, as shown in 
Fig. n^a), producing a two-dimensional (2D) potential 
perpendicular to the beam propagation, that varies spa- 
tially but is stable in time. A nice feature of the optical 
setup is that the path of the speckle beam is combined 
with the path of the imaging beam using a dichroic mir- 
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ror (A > 800 nm reflected, A < 800 nm transmitted). 
The imaging setup is used to detect both the BEC and 
the speckle pattern, and thus we can image in consec- 
utive photos the position of the trapped BEC and the 
precise reahzation of the random potential that the BEC 
experiences. The BEC is produced in a loffe-Pritchard 
magnetic trap, elongated perpendicular with respect to 
the speckle beam. The trapping frequencies are coz = 
27r X (8.74 ± 0.03) Hz axially and uj^ =27rx (85 ± 1) Hz 
radially, with the axis of the trap oriented horizontally. 
Our typical BECs are made of 3 x 10^ atoms in the 
hyperflne ground state \F = l;mi? = —1 >, with a peak 
density n 1.2 x 10^^ cm~^. 

A 3D representation of a typical speckle potential is 
displayed in Fig. ^b). In the right part Fig. ^ we show 
the Fourier transform of this potential, where the small- 
est length scale of the speckle is 10 /im. The average 
distance between neighboring speckles is approximately 
20 jam. The axial and radial size of our BEC is 110 /im 
and 11 /im respectively. The BEC probes 6 wells in the 
axial direction and only 1 well on average in the radial 
direction, and thus the cylindrical geometry enforces a 
quasi-one dimensionality on the speckle-condensate sys- 
tem. We deflne the speckle height Vs by taking twice 
the standard deviation of the speckle potential over a 
distance of 200 jam along the BEC axial direction. The 
speckle height is expressed throughout this Letter in units 
of frequency. 
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FIG. 2: (a) Density profile of the BEC after 28 ms of expan- 
sion from the combined magnetic and speckle potential for 
varying speckle potential intensities Vs (indicated in bottom 
part). To verify that the transfer to the speckle is adiabatic 
(b) shows the density profiles with (i) no speckle potential, (ii) 
when the speckle potential is abruptly switched off, and (iii) 
when it is adiabatically ramped off. (c) Dipole oscillations in 
the combined magnetic and speckle potential for varying Vs. 

An initial insight into the effect of the random poten- 
tial can be obtained by considering the BEC in expan- 
sion. Fig. I2Ia) shows a series of absorption images af- 



ter 28 ms of free expansion for different values of Vs. In 
each image the speckle intensity is increased adiabatically 
using a 200 ms exponential ramp with a time constant 
of 40 ms. After 50 ms in the combined magnetic and 
speckle potential, both potentials are abruptly switched 
off. For Vs ~50 Hz much lower than the chemical poten- 
tial of the BEC in the harmonic trap (~1 kHz), no change 
from the standard Thomas-Fermi (TF) proflle is seen. 
As the potential is increased to an intermediate power 
(~300 Hz) we see a strong modiflcation in the expanded 
density proflle. Such modiflcations can be a signature 
of the development of different phase domains across the 
BEC, similar to the stripes seen with highly elongated 
quasi-condensates [ij, , or when different momentum 
components are present due to the growth of instabilities 
17]. We note that the spacing of the stripes we observe, 
> 50 /im, can not be related to the characteristic length 
scale of the speckle pattern, 10-20 /im, the latter would 
give rise to an interferogram length scale of 6-12 jam after 
28 ms of expansion. Finally, when Vg >1 kHz, a value 
greater than the chemical potential, the system enters the 
tight-binding regime. Here in expansion only a broad 
gaussian proflle is seen. The tunnelling between wells 
is expected to be completely suppressed for such large 
spacing of adjacent sites (> 10 /im), and we are seeing 
the expansion of an array of randomly spaced sources. 
Unlike the situation of the expansion of a ID array of 
equally spaced sources with random phases where a pe- 
riodic interference pattern is still visible in a single run 
[T^ . no regular interference is expected to be seen with 
random spacing either for uniform or for random phases. 
To check that the ramp is adiabatic, we have ramped the 
speckle on (to our maximum speckle height of 1.7 kHz) 
and then off, using the exponential ramp described ear- 
lier. Fig. I2b)(i) shows the expanded BEC density pro- 
flle with no speckle, (ii) the expanded proflle when the 
speckle potential is abruptly switched off, and (iii) the 
recovered TF proflle when the speckle intensity is adia- 
batically decreased. This also clearly indicates that the 
density proflle in Fig. [2Ub)(ii) can not be attributed to 
heating of the atomic cloud. 

The transport properties in each of the above regimes 
are investigated by abruptly shifting the magnetic trap 
in the axial direction by 25 /im inducing dipole oscil- 
lations shown in Fig. El^c). Undamped oscillations are 
observed in the absence of the speckle. The dipole mo- 
tion for very weak Vs (~60 Hz) is slightly damped with an 
unchanged frequency, and in the intermediate 400 Hz po- 
tential, where phase fragmentation is seen in expansion, 
strong damping is observed. Once in the tight-binding 
regime, Vs >l kHz, the atomic cloud does not oscillate 
and remains localized in the deep, largely spaced, wells 
of the speckle potential on the side of the magnetic trap. 
This is analogous to the 'pinning' effect seen in magnetic 
microtraps, where the axial frequency is relaxed but the 
atoms do not propagate along the newly created waveg- 
uide remaining localized in the local minima due to the 
disorder in the microchip ^13.] . 
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FIG. 3: Quadrupole oscillations with increasing Vs for two 
different speckle realizations: (a) frequency shift (relative to 
the quadrupole frequency 13.72 zb 0.06 Hz measured in the 
pure harmonic trap), and (b) damping rate. The dashed line 
corresponds to \/^/2 Ud, calculated from the measured dipole 
frequency of 8.74 zb 0.03 Hz. The error (dotted lines) reflects 
the variation of the trap over several months. 



Precise information on the nature of BEG has been 
gained through studies of collective excitations 19] and 
it is therefore a natural starting point to gather infor- 
mation on a BEG in a random potential. We measure 
the frequency of both the dipole, and quadrupole, Vq, 
modes of a BEG confined in the harmonic magnetic trap 
combined with a weak speckle potential in the range of 30 
- 200 Hz corresponding to a small perturbation of the har- 
monic potential. The quadrupole mode is excited with 
a resonant modulation of the radial trapping frequency 
at 13.7 Hz, producing oscillations of 5% of the trapped 
BEG axial width, while the dipole mode is excited by an 
abrupt shift of the trap of 25 jam in the axial direction. 
In both cases the speckle potential remains stationary. 
The results for the frequency shift of the quadrupole os- 
cillations versus Vs for two different speckle realizations 
are shown in Fig. [31 with the corresponding damping of 
this mode. The shift is taken relative to the measured 
quadrupole frequency, 13.72 ± 0.06 Hz, in the pure har- 
monic trap. We see frequency shifts of up to 3%, and also 
observe that both the sign and amplitude of the shift de- 
pend on the exact realization of the speckle potential. 
We measure no change in the dipole frequency within 
our experimental resolution of 1%. These measurements 
are in contrast to what is seen with a regular standing- 
wave optical lattice, where both and Uq are rescaled 
with the single particle effective mass due to the optical 
lattice 20J. 

In the absence of the speckle potential, the quadrupole 
mode of an elongated BEG is given by Un = \J~^j2vd in 
the zero-temperature TF approximation ^| . In the pres- 
ence of an additional potential, this value can be modi- 
fied by both the deviation of the potential from a pure 
harmonic oscillator, and other effects such as changes in 
the condensate and superfluid fraction, modification of 
the interatomic interactions, or an increase in the ther- 
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FIG. 4: (a) Dipole and quadrupole frequency shifts for 150 
slices of the same speckle realization calculated from the sum 
rules approach, (b) Measured shift in the quadrupole fre- 
quency for 11 different realizations of the speckle at a fixed 
height of 80 =b 10 Hz. 



mal component. For example, in the collisional hydrody- 
namic limit, calculated for a gas above the condensation 
temperature, the frequency of the q uadrupole mode in a 
harmonic potential is Vq = ^JVl|hvd, while in the non- 
interacting thermal limit Vq = 2vd Actually, even 
without the speckle potential, our measured quadrupole 
frequency is shifted —0.7% with respect to the predicted 
^Jb/2h'd^ similar to what has been measured previously 
and attributed to the presence of a residual thermal com- 
ponent To check that the observed frequency shifts 
with the speckle present are not simply an effect of heat- 
ing, we have measured the quadrupole frequency in the 
presence of a large thermal fraction (50%) observing only 
a small shift of -0.4% with respect to the measurement 
with no discernable thermal fraction. 

To understand the effect of a shallow random potential 
on the dynamics of the BEG we use two theoretical ap- 
proaches. The Gross-Pitaevskii equation (GPE) is solved 
in the actual potential used in the experiment, and the 
numerical predictions are confirmed using a sum-rules 
approach. We investigate the dynamics of the BEG with 
the GPE in the combined speckle, Vopti and harmonic, 
V/io, potential 



dt 



2m 



opt 



Vho 



^(1) 



where ^ is the complex BEG order parameter, m is the 
atomic mass, and a is the scattering length. The GPE is 
solved after a sudden displacement Az of the harmonic 
trap for the dipole mode or after a sudden change in the 
axial trapping frequency for the quadrupole mode, for a 
variety of different speckle realizations. In the case of 
small amplitude oscillations {Az = Sjam and quadrupole 
oscillations of 5% of the BEG aspect ratio), the BEG os- 
cillates coherently with no appreciable damping for sim- 
ulations up to 600 ms. The frequencies are red or blue 
shifted depending on the particular speckle realization, 
and in some configurations a red shift is observed for 
one mode and a blue shift for the other. The result for 
the quadrupole mode shows an increasing frequency shift 
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with increasing Vg as qualitatively confirmed by the ex- 
perimental results in Fig. [SJi. We found that the pre- 
dicted frequency shift of the dipole mode reduces for 
larger amplitude oscillations since in this case the BEC 
probes more of the outer unperturbed harmonic oscilla- 
tor potential, and less of the central perturbed potential. 
Significant frequency shifts are observed (> 1%) only for 
small amplitude oscillations {Az < 10/im), but unfortu- 
nately we do not have sufficient experimental accuracy 
with the small amplitude oscillations to measure these 
shifts. The somewhat surprising behavior predicted by 
the GPE of disparate dipole and quadrupole frequency 
shifts can be explained by using the sum rules approach 
and considering Vopt as a perturbation. According to the 
sum rules, the frequencies of the low lying collective ex- 
citations can be estimated as h'^uj'^ = ms/mi where the 
moments rrii can be expressed as commutators involving 
the many-body hamiltonian H and a suitable excitation 
operator F In our case this operator can be cho- 

sen diS Fd ^ z and Fq ^ r'^ — az^ for the dipole and 
quadrupole modes respectively, a being a variational pa- 
rameter 24J. Treating the speckle potential Vo'pt as a 
small perturbation, and writing cc;^ = cc^q + A we get 

^ ^R'Kp.)o (2) 

m {z^)o 

where the second line is obtained assuming a strongly 
elongated BEC. These analytical results demonstrate 
how the dipole and quadrupole frequencies have differ- 
ent dependencies on the particular characteristics of the 
perturbing potential. We have verified that these predic- 
tions are in agreement with the full solution of the GPE 
in the cases considered. 

In Fig. 2fa) we show the predictions of the sum rules 
for 150 realizations of the speckle potential, with Vs ~ 
30 Hz. Each realization is obtained by taking succes- 
sive ID slices of 1 jam from the same 2D experimental 



speckle pattern. The average shift is close to zero for 
both excitation modes. We have performed a similar se- 
ries of measurements experimentally. Fig. ^b) shows 
the quadrupole frequency for different realizations of the 
speckle potential at Vs ~80 Hz. As predicted, we ob- 
serve both red and blue shifts, however within our limited 
statistics we see a bias towards the red. We have also be- 
gun an exploration of the dipole and quadrupole modes 
for larger amplitude oscillations where the solution of the 
GPE indicates that the superfluidity of the system may 
be compromised. The simulations show short wavelength 
modulation of the density distribution and damped os- 
cillations. Also in the experiment we found first evidence 
of increased quadrupole frequency shifts and strong per- 
turbations of the expanded density profile in this regime. 

In summary, we present the first results on a BEC in a 
controllable random speckle potential. We see the effect 
of a weak random potential indicated by stripes in the 
expanded density profile of the BEC, damped dipole os- 
cillations, and frequency shifts in the quadrupole mode 
that are not correlated to the measured dipole frequency. 
Simulations solving the GPE and an analytical sum-rules 
approach successfully describe these shifts in the pertu- 
bative regime. In the limit of a strong random potential, 
we observe a broad gaussian density profile after expan- 
sion, indicating that the atoms remain localized in the 
individual wells of the speckle potential. This opens the 
way to new investigations of localization phenomena in 
quantum gases, and the superfluid behavior with varying 
roughness of a surface potential. 
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